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1.  Introduction 


In  order  to  make  the  photo-acoustic  spectrometer  used  for  the  detection  of  energetics  more 
efficient,  it  is  important  to  understand  the  acoustic  resonance  conditions  in  the  photo-acoustic 
cell.  At  first  glance  this  may  seem  simple,  but  as  it  turns  out  the  interaction  of  the  main  tube  with 
side  branches  is  quite  complicated.  Assuming  that  a  hole  is  made  exactly  in  the  middle  of  the 
main  tube,  the  length  and  diameter  of  the  side  branch  will  determine  the  resonant  frequency.  This 
can  be  anywhere  between  the  resonance  of  the  main  tube  and  twice  that  frequency.  Starting  as 
early  as  1922  and  continuing  through  the  present  day,  researchers,  such  as  Mason  (i)  at  Bell 
Labs  have  studied  the  acoustics  of  tubes  and  tried  to  develop  techniques  that  would  predict  the 
resonant  or  filtering  characteristics  of  tubes  with  side  branches. 

As  long  ago  as  the  Golden  Age  of  Greece,  man  has  been  able  to  place  holes  in  the  correct 
position  along  a  hollow  tube  to  make  music.  Today  instrument  designers  have  developed  a 
mathematical  technique  to  determine  the  location  and  shape  of  the  holes  along  the  tube. 

A  more  modem  approach,  which  is  an  extension  of  Mason’s  use  of  electrical  network  theory 
applied  to  acoustics,  is  the  use  of  electrical  transmission  line  theory  to  describe  the  tube  and  the 
side  branch  as  matrices.  The  resultant  matrix  is  then  applied  to  the  pressure  and  mass  flow  at  any 
point  along  the  tube  in  order  to  calculate  the  pressure  and  mass  flow  at  any  other  point.  Here, 
pressure  and  mass  flow  are  synonymous  with  voltage  and  current. 

In  this  report,  I  will  discuss  both  musical  instrument  design  and  matrix  theory  that  is  applied  to 
acoustic  tubes  with  side  branches.  I  will  also  discuss  a  series  of  experiments  performed  on  tubes 
with  side  branch  holes  of  various  sizes  and  shapes.  I  will  then  apply  the  theory  to  our  photo¬ 
acoustic  spectrometer  to  show  that  the  theory  predicts  the  resonant  frequency  within  a  few 
percent.  In  closing,  I  will  discuss  several  sources  of  error  in  these  methods. 


2.  Acoustic  Waves  in  Cylindrical  Tubes 


2.1  Plane  Waves 

Acoustic  waves  that  are  confined  to  a  tube  with  rigid  walls  can  be  considered  to  be  plane  waves, 
assuming  the  diameter  of  the  tube  is  small  compared  to  the  wave  length  of  the  acoustic  wave. 
Higher  order  modes  are  only  possible  when  the  diameter  of  the  tube  is  greater  than  two-thirds  of 
the  wave  length  (2).  Also,  the  plane  wave  approximation  only  holds  for  frequencies  where  the 
effects  of  viscose  friction  are  negligible.  At  higher  frequencies,  the  effect  of  viscosity  is  confined 

(2v\1!2 

to  a  thin  layer  that  is  close  to  the  wall  ( 8a  =  ( —  1  ),  where  v=q/p  is  the  kinetic  viscosity.  In 
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order  to  make  the  plane  wave  approximation  valid,  one  should  have  a  boundary  layer  that  is  thin: 
8 

«  1  (where  d  is  the  diameter  of  the  tube).  This  leads  to  the  following  relationship  (2): 

«  f  <—.  The  kinetic  viscosity  of  air  is  1.5  x  10_5m2/s,  so  for  a  representative  tube  with  1 

J  L  LA*  LA* 

mm  diameter,  the  plane  wave  approximation  is  valid  for  frequencies  9.5  Hz<f<170  kHz. 

When  acoustic  waves  propagate  (alternating  compression  and  expansion  of  the  medium),  the 
temperature  gradients  between  these  compressed  and  expanded  volumes  are  very  small.  Because 
of  this,  very  little  heat  flows  from  one  region  to  another,  and  the  process  can  be  considered 
adiabatic.  In  an  adiabatic  process,  the  pressure  in  the  tube  can  be  represented  as  a  function  of 
density,  P=P(p).  Now  from  the  ideal  gas  law  PV=nRT,  and  for  an  adiabatic  process,  this  leads  to 
P=Cpy,  where  C  is  a  constant  and  y  is  the  ratio  of  the  specific  heats.  If  we  now  represent  the 
function  P(p)  as  a  Taylor  series  expansion  around  the  point  po  and  discard  higher  order  terms  we 
have 


p-p«=Otp-p^ 

(1) 

In  the  limit,  this  can  be  expressed  as 

dP=O0d»’ 

(2) 

where  ( — )  is  the  slope  of  the  adiabatic  plot  of  pressure  vs.  density  at  the  point  (Pofio)- 

\dp/  q 

Acoustic  waves  are  longitudinal  waves.  This  means  that  particles  move  back  and  forth  in  the 
direction  of  propagation,  creating  alternating  volumes  of  compression  and  expansion.  The 
restoring  force  that  propagates  the  wave  along  is  simply  the  elastic  opposition  that  arises  when 
the  fluid  is  compressed. 

Figure  1  below  will  be  used  to  derive  an  acoustic  plane  wave  equation.  The  variables  to  be  used 
in  this  derivation  are  listed  and  defined  below  the  figure.  This  derivation  follows  the  discussion 
given  by  Kinsler  and  Frey  ( 3 ). 
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Figure  1.  Diagram  of  a  slice  down  the  axis  of  the  tube. 
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The  variables  to  be  used  in  this  derivation  are  listed  and  defined  below: 


x=  equilibrium  coordinate  of  a  thin  sheet  of  molecules 
=  displacement  of  the  sheet  along  the  axis  of  propagation 
u  =  velocity  of  the  sheet w  =  — 

p  =  instantaneous  density  at  any  point 

po  =  constant  equilibrium  density  of  the  medium 

5  =  condensation  at  any  point  defined  below 

s  =  Pz£lorp=p0(i+s)'  (3) 

Po 

P  =instantaneous  pressure  at  any  point 
Po  =  constant  equilibrium  pressure 

p  =  excess  pressure  or  acoustic  pressure  at  any  point  defined  as  p  =  P-Po 
c  =  velocity  of  wave  propagation 

First,  we  must  address  the  conservation  of  mass.  The  mass  contained  between  x  and  x+dx  can  be 
written  as  poSdx,  where  Sdx  is  the  volume  element.  After  an  acoustic  wave  passes,  the  particles  at 
x  will  have  moved  to  x+^,  and  the  particles  at  x+dx  will  have  moved  to  x+dx+(l+dc),  so  that  the 
new  volume  is  Sdx(\+-^).  The  mass  within  this  volume  must  be  the  same  as  the  original  mass 
due  to  the  conservation  of  mass,  which  leads  to  the  following: 

pSdx(\+j-x)  =  poSdx.  (4) 

Now,  if  we  replace  p  with  po(\+s)  from  equation  3,  equation  4  becomes 

(l+s)(l+g)=l.  (5) 


Because  both  the  density  and  molecular  displacements  are  very  small  (for  even  load  noises 
neither  will  exceed  10  4),  we  can  neglect  the  product  of  s^,  and  equation  5  becomes 


s  = 


dj_ 

dx 


(6) 


For  acoustic  waves,  the  incremental  pressure  changes  are  very  small,  so  that  dP  can  be  replaced 
by  p  the  acoustic  pressure,  and  the  small  change  in  density  dp  can  be  replaced  by  spo  (see 
equation  3).  Now,  using  these  in  equation  2  leads  to 


P  = 


(7) 
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We  can  also  let 


so  that  we  have 


p  =  poc2s. 

Finally,  by  replacing  s  with  —  ^  from  equation  6,  we  have 

2  a? 

P  =  ~Poc Zyx- 


(8) 

(9) 


(10) 


When  a  fluid  is  deformed,  as  we  have  described,  the  pressure  acting  on  the  two  faces  of  the 
volume  element  Sdx  is  slightly  different,  producing  a  net  force  in  the  positive  x  direction. 

Because  the  external  force  acting  on  each  face  is  equal  to  the  pressure  times  the  area,  the  net 
force  can  be  written  as 

dFx=[p-  (v+£dx)]s=-d£dxS.  (11) 


Now,  if  we  set  this  net  force  equal  to  the  mass  of  the  volume  element  times  its  acceleration,  we 
are  left  with 


dp  _  a2f 
~dx~  P°  at2 ' 


(12) 


Differentiating  equation  10  with  respect  to  x  and  replacing  ^  in  equation  12,  we  are  left  with 


d2l=c2d_ll 

dt2  dx2‘ 


(13) 


This  is  the  wave  equation  in  displacement,  and  similar  equations  apply  to  pressure,  velocity,  and 
condensation. 


A  note  here  should  be  added  in  regard  to  the  wave  velocity  in  a  tube.  In  air,  the  velocity  is  a 
function  of  temperature  alone  and  is  given  by 

c  =  331.3  +  .606T(°C)  j.  (14) 

At  20  °C,  this  leads  to  a  velocity,  c  =  343.6 -j.  This  velocity  is  no  longer  valid  inside  a  tube  with 

rigid  walls  due  to  the  thermal  exchange  with  the  walls  and  viscous  friction  at  the  walls.  The  first 
to  study  the  velocity  and  attenuation  of  sound  waves  in  tubes  was  Helmholtz  ( 4 )  in  1863.  Later, 
Kirchhoff  (5)  and  Rayleigh  (6)  extended  this  theory.  Zwikker  and  Kosten  (7)  offered  the  most 
widely  accepted  equation  for  the  propagation  constant  derived  from  Kirchhoff  s  work: 
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(15) 


I  = 


lMi3/lS)  I 

V  M^S)  V 


~!  +  (-) —  1) 


Where  Jn  are  the  Bessel  functions  of  order  n;  S  =  r ^j—  is  the  shear  wave  number;  r  is  the  inner 

radius  of  the  tube;  co  is  the  radian  frequency;  the  parameter  \[oS  is  the  ratio  of  the  radius  r  to  the 

/2v\V2 

thermal  boundary  layer  defined  earlier  as  (  8a  =  I  —  1  ),  where  v=r|/p  the  kinetic  viscosity;  and 

Q 

Y  =  -jr  is  the  ratio  of  the  specific  heats.  Now,  the  important  thing  to  take  away  here  is  that  the 
phase  velocity  is  given  by 


v„  = - 

p  lm(r) 


(16) 


Tijdeman  (8)  did  a  numerical  analysis  of  this  equation  and  published  the  result  in  tables  based  on 
S  and  k\  the  reduced  frequency  given  by  k  =  — .  For  our  representative  cell,  with  1  mm  diameter 

and  10  mm  length,  S= 44  and  Ar=0. 1 69 1 .  From  the  tables,  the  Im(f)  =  1.0261  and  the  velocity  in 
the  tube  is  335  m/s. 


2.2  Higher  Order  Modes  (9, 10, 11) 

Figure  2  depicts  a  section  of  a  cylindrical  tube  with  radius  r=a.  It  is  assumed  that  the  tube  walls 
are  perfectly  rigid  and  smooth.  Wave  propagation  is  along  the  positive  z  axis  into  the  report.  We 
find  that  the  acoustic  pressure  in  the  tube  is  governed  by  the  wave  equation: 

1  d 2 

V2p(r,(p,z,t)--—p(r,(p,z,t )  =  0.  (17) 


Figure  2.  Geometry  of  a  tube  in  cylindrical  coordinates. 

We  can  use  the  technique  of  the  separation  of  variables  to  solve  this  differential  equation  by 
allowing 

p(r,(p,z,t )  =  R(r)&((p)Z(z)T(t).  (18) 
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The  Laplace  operator  in  cylindrical  coordinates  is  defined  as 


a2  i_a_  ,  _d_ 

Qr2  r  Qr  r2  Q(p2  Qz2 


)■ 


(19) 


If  we  now  insert  equation  1 8  into  equation  1 7  and  apply  the  Laplace  operator,  we  are  left  with 
the  following  after  dividing  through  by  p(r,  cp,  z,  t): 

R"(r)  lR'(r)  1  <?"(<?)  j  Z"(z)  1  T"(t)  _ 

7?(r)  r  fl(r)  r2  <P((p)  Z(z)  c2  7(t)  *  ^ 


In  order  to  proceed  with  the  solution  of  equation  20,  move  the  term  exclusive  in  t  to  the  right- 
hand  side  of  equation  20.  Noticing  that  the  right-hand  side  is  a  function  of  t  alone,  it  must, 
therefore,  be  equal  to  a  constant.  Allowing  that  constant  to  be  equal  to  —  k2,  where  k  =  — ;  we  are 

led  to  the  following  equation: 

T"(t)  +  co2T(t)  =  0.  (21) 


This  equation  is  that  of  a  harmonic  oscillator,  and  its  solution  is  of  the  form 

T(t)  =  T0e±j(Ot. 


(22) 


Now,  we  address  the  term  in  z.  We  can  move  the  term  in  z  to  the  right  side  of  equation  20  so  that 
now  the  right  side  of  equation  20  is  a  term  in  z  alone.  The  right  side  again  must  be  equal  to  a 
constant  as  follows: 

-fc2-S=-^-  <23> 


This  can  be  rearranged  as 


Z"(z)  +  ( k 2  -  a2)Z(z )  =  0. 


(24) 


The  solution  of  this  equation  will  be  of  the  form 

Z(z)  =  Z0e±jPz.  (25) 

Here,  /?2  =  k2  —  a2.  It  should  be  noted  that  when  k2  >  a2,  the  solution  is  oscillatory;  when 
k2  <  a2,  the  solution  is  exponentially  growing  or  decaying;  and  when  k2  =  a2,  the  solution  is 
constant.  By  convention,  we  will  choose  k2  >  a2,  and  if  we  choose  propagation  in  the  positive  z 
direction,  we  can  drop  the  +  sign  from  the  exponent. 


Equation  20  can  now  be  written  in  the  form 


r2  d2R  ,  r  dR  ,  99 

- 7  + - 1- 

R  dr2  R  dr 


-1  d20 
0  dcp 25 


(26) 


and  because  we  have  two  completely  independent  equations  in  different  variables  equal  to  each 
other,  they  must  both  be  equal  to  the  same  constant  so  that 
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-1  d2^  _  2 

<P  d<f 2  _  <T<iP’ 

or  0  +  =  (27> 

and  r2  +  r  ^  +  (jV  -  <r$)R  =  0.  (28) 

The  solution  to  equation  27  is 

0((p)  =  a0sin(o(p(p)  +  b0cos(a<p(p),  (29) 

and  0((p)  must  be  periodic  in  <p  or  O (<p)  =  <P((p±m2n),  so  that  o(p  must  be  an  integer  n=l,2,3. . 

O(cp)  =  a0sin(mp)  +  b0cos(mp).  (30) 

Now,  equation  28  is  Bessel’s  differential  equation.  Solutions  to  this  equation  take  the  form  of 

R(r)  =  CJn(ocr)  +  DYn(ocr),  (31) 

where  J„  and  Yn  are  the  Bessel  functions  of  the  first  and  second  kind  in  order  n,  respectively. 
There  are  two  boundary  conditions  that  must  be  met.  First,  the  function  p(r,(p,z,t)  must  be 
bounded  everywhere  within  the  tube.  Because  the  Bessel  function  of  the  second  kind  approaches 
—oo  as  r  -*  0,  we  must  set  D  =  0.  The  second  boundary  condition  is  as  follows.  Because  the 
wall  of  the  tube  is  perfectly  rigid,  there  can  be  no  net  flow  in  the  direction  of  r  at  the  wall  r  =  a: 

^  =  0  at  r  =  a.  (32) 

dr  v  y 

We  can  now  state  the  complete  solution  of  the  wave  equation  for  p(r,tp,z,t): 

p(r,(p,x,t)  =  (Asin  (ncp)  +  Bcos(n<p))Jn(acr)exp[j(-fiz+GOt)] ,  (33) 

where  the  constants  C,  Zo,  and  To  have  been  incorporated  into  A  and  B.  The  boundary  condition 
at  r  =  a  leads  to  J’n(aca)=0.  After  looking  up  the  roots  of  the  first  derivative  ]'n  ( aca  )  =  0,  we 
find  that  cca  =  1.841 18  for  the  first  propagating  none  0,0  mode  of  the  lowest  order.  Because  ac  = 
2n/Xc  and  A  =  c/fc,  where  c  is  the  velocity  of  sound,  and  fc  is  the  cutoff  frequency: 

1.841 18  =  2na/Xc  =  2naf/c, 

or  fc  =  1.841 18c/Ara.  (34) 

If  we  apply  this  to  a  representative  tube  with  1  mm  diameter,  and  we  assume  an  approximate 
velocity  of  340  m/s  fc  =200  kHz,  this  indicates  that  any  wave  of  frequency  less  than  200  kHz 
will  propagate  as  a  plane  wave. 


7 


3.  Transmission  Line  Theory  Applied  to  Acoustics 


It  turns  out  that  a  duct  with  an  acoustic  wave  propagating  within  can  be  modeled  as  if  it  was  an 
electrical  transmission  line.  In  all  that  will  follow,  this  analogy  will  be  useful.  The  analogous 
terms  are  pressure  with  voltage  and  mass  flow  with  current.  Acoustic  impedance  will  be  defined 
as  pressure/mass  flow  (Z  =  ^  =  ^),  and  admittance  is  defined  as  K  The  density  of  the  gas  is  po, 

c  is  the  speed  of  sound  in  the  medium,  and  S  is  the  cross-sectional  area  of  the  tube.  Expanding  on 
the  transmission  line  analogy,  distributed  inertance  (M)  and  distributed  compliance  (C)  are 
associated  with  inductance  and  capacitance,  respectively.  As  in  a  transmission  line,  the 
characteristic  impedance  is  given  by  (12) 


where  M  =  —  and  C  =  ———  . 

S  p0c2 

Now  let’s  begin  the  analysis  of  a  transmission  line  using  figure  3  as  our  model. 


Po/S 

S/p0c2 


PqC 

S 


(35) 


A  transmission  line  is  different  than  an  ordinary  electrical  circuit,  because  its  components  are 
distributed  along  the  line  and  not  lumped  into  resistors  and  capacitors.  It  turns  out  that  if  the  load 
impedance  is  not  the  same  as  the  characteristic  impedance,  there  will  be  a  reflected  wave 
traveling  back  toward  the  source  due  to  the  conservation  of  energy  at  the  load.  In  this  case,  that 
direction  is  left.  It  also  turns  out  that  the  impedance  changes  along  the  line.  This  is  due  to  the 
load  impedance  having  a  differing  effect  at  different  positions  along  the  line.  This  is  analogous  to 
the  way  the  load  on  the  secondary  of  a  transformer  is  seen  by  the  primary  due  to  mutual 
inductance  or  coupling  between  the  two  sides: 


Z(z)  = 


nz) 

I(z)' 


(36) 


To  develop  an  equation  for  the  impedance  along  the  line  and  the  input  impedance  seen  by  the 
source,  we  have  to  look  at  the  superposition  of  the  forward  and  reverse  voltage  and  current: 
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V{z)  =  V+e~jkz  +  V_ejkz  ,  (37) 

7(z)  =  Y0(V+e~Jkz  -  V_ejkz).  (38) 

Here,  Yo  is  the  admittance  into  the  transmission  line  (Yo=l/Z),  V+  is  the  wave  traveling  right,  and 
V.  is  that  traveling  left.  Figure  4  shows  a  segment  of  the  transmission  line.  If  we  apply 
Kirchhoff  s  voltage  law  around  this  loop,  we  have 

V  -  ja)Ll{z)dz  -  (V  +  dV)  =  0,  (39) 

from  which  we  have  ^  =  — /o)L/(z).  (40) 

Using  the  fact  that  k  =  oW LC ,  (41) 


and  applying  equation  40  to  equation  37,  we  are  left  with  equation  38. 


Now,  to  develop  an  expression  for  the  impedance  as  a  function  of  position  on  the  line,  we  need  a 
couple  of  items.  First,  Z[  =  V(z=l)/I(z=l).  Then  we  need  to  define  a  reflection  coefficient  f i=V. 
ZV+.  If  we  apply  these  to  equation  37  divided  by  equation  38  and  let  the  load  be  at  z=0  we  have 


y++y_  _  v++rLv+  v  _  i+rL  v  _  v 
i+-i-  v+-rLv+  u  l -rL  u  L 


(42) 


Solving  this  for  fz,  we  have 

j-i  _ Zl—Zq 

L  z0+zL 

Zfz)  =  ^  zz  v+e~jkz+rLV+ejkz 
('ZJ  Kz)  -L(V+e-ikz-rLV+eikzy 


(43) 


7(ys  7  (Z0+ZL)e-tfcz+(ZL-Z0)effcz 
^  }  0  (Z0+ZL)e-ffcz-(ZL-Z0)efkz’ 
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Z(z)  =  Z0 


ZL  cos(kz)-jZ0  sin(kz) 


Z0  cos(kz)-jZL  sin(kz)' 

Letting  the  position  of  the  source  be  z=-/,  we  are  left  with  the  important  result 

7  _  7  ZL+jZ o  tan(kt) 

0  Z0+ tan(ki)' 

This  equation  will  be  used  repeatedly,  as  we  develop  the  equations  for  resonance  in  our  tube. 


(44) 


v2 

A 

B 

Vi 

h 

C 

D 

h 

Now  we  apply  a  very  powerful  matrix  representation  to  our  transmission  line  problem.  We  want 
to  develop  a  technique  that  will  allow  us  to  transform  a  voltage  and  current  at  any  point  along  the 
line  to  a  voltage  and  current  at  a  different  point  on  the  line: 


(45) 


Let’s  be  reminded  of  equations  37  and  38,  and  we  will  let  position  1  be  at  z=0  (the  load)  and 
position  2  be  at  z=-l  (the  source)  as  we  have  discussed  before: 

V(z)  =  V+(e~jkz  +  rLejkz ),  (46) 

/ (z)  =  Y0 7+  (e ~jkz  -  rLejkz).  (47) 

Now  let’s  let  the  load  be  an  open  so  that  rL  =  +1  and  V(z=0)=+2V+  applying  these  to  equations 
45,  46,  and  47  we  have 


These  lead  to 


V+(ejkl  +  e~jkl )  =  2AV+, 

(48) 

Y±(eW  -  e-JM)  =  2 CV+. 
z0 

(49) 

A  =  cos  kl, 

(50) 

C=  —  sin  kl. 

Zo 

(51) 

Conversely,  if  we  let  the  load  be  a  short  so  that  f" L=  -1  and  l(z=0)=—  V+  and  V(z=0)  =  0  again, 

zo 

using  equations  45,  46,  and  47  we  have 


V(z  =  - 1 )  =  V+(eJkl  -  e~jkl )  =  —V+B, 

Zq 

i(z  =  -0  =  vf{e>kl  +  e~>kl)  =  f  V+D, 

z0  z0 


(52) 


leading  to 


B  =  jZ0  sin  kl. 


(53) 
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D  =  cos  kl. 


(54) 


So,  we  now  have  the  transfer  matrix  from  point  1  to  point  2,  where  /  represents  the  length  of  that 
section  of  the  line,  and  Zq  is  the  characteristic  impedance  in  that  section. 

We  can  extend  this  technique  to  handle  a  line  with  any  number  of  sections  in  the  following 
manner.  We  define  the  transfer  matrix  as  follows: 


Ti  = 


cos  kill 
—  sin  k^i 

zoi 


jZoi  sin  kili 
cos  k^i 


(55) 


and  we  use  this  to  describe  a  line  with  N  sections  as  follows: 

I  Vi, 


-  tn  ■  tn_1  •■■T1 


I, 


(56) 


This  will  become  useful  later  when  we  describe  acoustic  tubes  with  side  branches  and  tubes  that 
change  diameter. 


4.  Musical  Instrument  Design 


A  flute  is  essentially  a  hollow  tube  with  a  way  at  one  end  to  modulate  the  pressure  of  the  air 
column  within  the  tube.  Along  the  length  of  the  tube  are  holes  that  essentially  change  the  length 
of  the  tube,  and,  thereby,  change  the  frequency  of  resonance  in  the  tube.  We  will  start  with  a  tube 
with  no  holes  as  shown  in  figure  5.  At  resonance  in  the  fundamental  mode  a  differential  pressure 
standing  wave  is  established  in  the  tube  that  has  the  greatest  magnitude  in  the  center  of  the  tube 
and  falls  to  zero  at  the  ends.  Higher  order  plane  wave  modes  do  propagate,  and  these  will  always 
be  an  integral  number  of  half-wave  lengths.  Note  that  the  effective  length  of  the  bore  is  always 
longer  than  the  actual  physical  dimension  of  the  tube.  This  is  because  the  pressure  wave 
penetrates  a  short  distance  beyond  the  end  of  the  tube  before  partially  reflecting  back  into  the 
tube.  I  will  state  here  without  discussion  that  this  correction  is  0.3  Dbore  (3,  13,  14).  So  that  the 
effective  length  of  the  tube  is  Leffbore  =  Lbore  +  0.6DbOre  and  the  wave  lengths  of  the  resonant 
modes  are  given  by  nX  =  2Leffbore,  where  n  is  an  integer  (1,2,3,. . .). 


11 


Figure  5.  Tube  with  no  holes. 

The  addition  of  a  side  branch  effectively  changes  the  length  of  the  tube,  so  in  a  flute  when  one  of 
the  tone  holes  is  opened  the  sound  pitch  becomes  higher  because  the  effective  length  of  the  tube 
has  been  shortened.  In  what  follows,  I  will  discuss  the  effect  of  a  single  tone  hole,  or  side  branch, 
on  the  resonant  frequency  of  the  tube.  Figure  6  will  be  used  in  the  discussion  that  follows,  which 
is  based  on  that  of  Nederveen  (14). 


Figure  6.  Model  of  tube  with  singe  side  branch. 


As  discussed  in  section  3,  an  acoustic  tube  can  be  modeled  as  if  it  is  an  electrical  transmission 
line.  The  characteristic  impedance  (Z0)  is  equal  to  p/U,  the  differential  pressure  (or  acoustic 
pressure)  divided  by  the  mass  flow  associated  with  the  particle  movement  as  the  wave  passes. 
Now  applying  Kirchhoff  s  current  law  at  the  junction,  where  the  main  tube  meets  the  side 
branch,  we  see  that  the  current  into  the  junction  must  equal  the  current  out  of  the  junction. 
Because  the  pressure  for  all  three  branches  must  be  equal,  this  leads  to 


i 

ZL 


(57) 


or  the  admittance  looking  back  down  the  bore  from  the  side  branch  must  equal  the  sum  of  the 
two  admittances  at  the  junction.  Now,  if  we  return  to  equation  44,  we  can  describe  the  input 
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impedance  to  any  one  of  these  sections  by  first  setting  Zr  =  0,  because  all  of  our  tube  ends  are 
open,  and  acoustic  pressure  there  must  go  toward  zero.  We  are  left  with  the  following: 

Zjn  =  jZq  tan  kl.  (58) 

We  are  going  to  use  this  expression  many  times  with  modifications  for  the  length  of  the  tube 
being  described  and  the  characteristic  impedance  in  that  tube. 

Now  the  concept  of  a  replacement  tube  (14)  needs  to  be  introduced.  The  replacement  tube  is  a 
fictitious  tube  that  will  represent  the  resonant  tube  after  the  correction  for  the  side  branch  is 
introduced.  The  length  of  this  fictitious  tube  is  Leff  =  Le/-^bore  +  xb.  Here,  A/,  will  represent  the 

corrected  length  of  the  side  branch  te  adjusted  to  represent  the  side  branch  as  if  it  were  in  the 
main  tube  by  correcting  for  the  different  areas: 


Referring  to  the  expression  given  at  the  beginning  of  section  3  for  the  characteristic  impedance 
of  any  tube  and  equation  58,  we  can  express  the  general  admittance  into  a  tube  as  follows: 

Y=i=-(it^cot(kz+*y  (60) 

Here,  S  is  the  cross-sectional  area  of  the  tube,  and  t//  is  a  phase  introduced  to  handle  any  phase 
adjustment  needed.  Now,  in  order  to  describe  the  admittance  of  the  side  branch,  we  first  see  that 
at  z=0,  Y  ->  oo,  and  this  implies  that  \\i  =0,  so  that  at  z=-te: 

Ysb  =  (jfz)  cot  kte.  (61) 

\JPoCZ 

The  admittance  looking  down  the  tube  toward  the  load  end  is  found  in  a  similar  manner: 

Yr  =  cot/c^^I£  .  (62) 

Now,  to  develop  an  expression  for  the  admittance  looking  toward  the  source  from  the  side 
branch,  we  first  notice  that  at  resonance,  the  impedance  at  the  input  goes  to  zero,  so  that  Y  ->  oo. 
Allowing  z=0  at  the  source,  this  leads  to  \\i  =  -mn,  where  m  is  an  integer.  Now,  again  looking  at 
the  resonance  condition  for  the  replacement  tube,  we  have  kLejf  =  nm.  This  leads  to  the  following 

atz=_^effbor1. 

2 

Yl  =  ~  (f^)  cot (k  -  kLe„).  (63) 

Now,  we  can  develop  a  first  approximation  for  Xb,  which  will  be  used  in  describing  the 
experimental  results  to  follow.  Again  from  equation  57  we  have 

Yl  =  Ysb  +  Yr.  (64) 
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We  now  utilize  the  series  expansion  of  the  cot(x)  =  ^  because  for  our  approximation  we  drop 

all  the  terms  after  the  first.  Inserting  equations  61,  62,  and  63  into  64,  we  are  left  with  equation 
57  after  simplification: 


1 

Leff  bore 

2 


+  L 


£/^I£_Le// 


+ 


£ bare 
6  Ssb 


=  0. 


(65) 


Now,  we  can  let  h  represent  the  length  of  the  side  branch  as  it  would  appear  in  the  main  tube 
having  the  same  admittance  of  the  actual  side  branch: 

(66) 

$sb 


The  length  of  the  replacement  tube  used  to  describe  the  resonant  tube  of  length  Z^will  then 
given  by 


t  _  Leffbore  ,  1 

Leff  ~  2  h  Ab- 


(67) 


5.  Mason’s  Approximate  Network  Acoustic  Filter 


It  seems  that  Stewart  (15,  16)  was  the  first  to  describe  an  acoustic  filter  by  using  electrical 
network  theory  in  1922.  Mason  ( 1 )  showed  in  1927  that  this  theory  was  good  only  as  a  first 
approximation,  and  that  some  additions  had  to  be  made.  Mason’s  circuit  is  shown  in  figure  7. 


Here,  Zshunt  represents  the  impedance  of  the  side  branch;  it  could  contain  several  components.  Zq 
is  the  characteristic  impedance  of  the  transmission  line  leading  up  to  and  away  from  the  side 
branch.  In  Mason’s  case,  these  were  lumped  components.  This  circuit  is  very  similar  to  that 
which  is  still  described  today  by  modem  researchers.  For  this  reason,  I  am  including  it  here  for 
historical  perspective.  Later  in  the  analysis  section,  I  will  use  one  of  Mason’s  structures,  shown 
in  figure  8,  number  3,  and  his  equations  for  equivalent  length  and  area  to  calculate  the  resonant 
frequency  of  one  of  our  photo-acoustic  cells.  Mason  suggests  that  the  shunt  impedance  in  figure 
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7  should  contain  three  elements  in  a  T  or  n  configuration.  This  is  exactly  what  Keefe  ( 1 7) 
suggests,  and  is  in  fact  the  circuit  that  I  have  chosen  to  use  in  my  analysis. 


Figure  8.  Mason’s  equivalent  for  side  branch  representation. 


6.  Experimental  Results  and  Analysis 


Two  experiments  were  performed.  The  first  was  initiated,  because  we  were  trying  to  understand 
the  resonant  frequency  of  the  photo-acoustic  cell  in  our  photo-acoustic  spectrometer.  The 
performance  of  the  cell  during  regular  operation  always  seemed  to  have  a  resonance  above  the 
calculated  value.  There  were  two  reasons  for  this.  First,  it  is  very  difficult  to  measure  the  speed 
of  sound  in  the  cell.  In  section  2,  a  method  of  calculation  was  offered  that  yields  a  value  of 
335  m/s  in  the  tube.  This  is  the  value  I  will  use  in  the  calculations  done  for  our  1  x  10  mm  cell. 
Second,  we  were  not  accounting  for  the  length  correction  due  to  the  side  branch. 

The  first  experiment  was  done  to  simply  replace  our  cell  with  a  larger  version  made  of  PVC  pipe. 
This  tube  shown  in  figure  9  is  basically  a  scale  model  of  our  smaller  photo-acoustic  cell.  The 
photo-acoustic  spectrometer  cell  is  10  mm  in  length,  1  mm  in  diameter.  It  has  a  side  branch  that 
is  0.5  mm  in  diameter  and  1  mm  in  length.  The  model  cell  is  303  mm  long  with  diameter  of 
27  mm.  The  side  branch  is  20  mm  long  with  a  diameter  of  1 5  mm. 
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Figure  9.  Picture  of  experimental  setup  for  tube  with  a  single  side  branch. 
The  method  of  analysis  is  as  follows: 


Pin 

cos  kL  jZ0  sin  kL 

1  0 

1 

cos  kL  jZ0  sin  kL 

Pout 

Uin 

—  sin  kL  cos  kL 

z0 

—  1 

Zs 

—  sin  kL  cos  kL 

z0 

Uout 

These  matrices  represent  the  situation  shown  in  figure  10. 


-i  I  J, 

D»c,«  'X 

L  X  L  

\  /j\  / 

(68) 


Figure  10.  Schematic  of  the  tube  shown  in  figure  9. 
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In  equation  68,  Z$  is  the  shunt  impedance  that  represents  the  side  branch.  As  given  in  equation  58 


zs=jzo(^rs)  tan  fete . 


(69) 


The  term  ^p^orej  appears  because  the  characteristic  impedance  in  the  side  branch  is  inversely 

proportional  to  the  area  of  the  side  branch,  and  Zo  is  inversely  proportional  to  the  area  of  the 
main  tube.  After  calculating  the  total  matrix  for  equation  68  and  setting  Zin  =  — ,  we  are  left 

w  Uin 

with 


7.  — 

^ in 


jZ0( sin  2kL(Db°re}  tan kte+ (sin  2/cL)2 
V  V  Dsb  / 


(Dbore\ 
\  Dsb  ) 


(70) 


Now,  resonance  occurs  when  the  imaginary  part  of  Zin  is  equal  to  zero,  so  that 


sin2feL 


fPfcoreV 
V  Dsb  ) 


tan  ktP 


+  (sin  kL)2  =  0. 


(71) 


This  is  the  resonance  condition.  Note  the  term  in  the  square  brackets.  This  is  the  input  impedance 
to  the  side  branch.  Later,  when  we  analyze  side  branch  structures  that  are  different  than  the  one 
in  figure  10,  we  can  simply  calculate  the  new  input  impedance  to  the  new  structure  and  substitute 
it  in  the  square  bracket  in  equation  7 1 . 


Now  when  we  insert  the  dimensions  of  the  PVC  model  tube  into  equation  71  and  do  a  numerical 
analysis,  we  find  that  the  calculated  resonance  frequency  is  653  Hz.  In  the  test  done  on  the  tube 
in  figure  10,  the  measured  resonance  was  677  Hz.  This  is  an  error  of  3.7%.  The  same  analysis 
was  used  with  substitutions  for  the  dimensions  of  the  1  x  10  mm  cell.  The  calculated  resonance 
was  at  18,362  Hz.  The  measured  value  of  the  resonance  tested  with  vinyl  acetate  as  the  gas  in  the 
cell  was  18,500  Hz.  This  is  an  error  of  less  than  1%. 


A  second  experiment  was  performed  on  the  PYC  tube  in  figure  11.  We  started  this  test  with  a 
tube  that  had  no  hole  in  the  side.  We  tested  the  resonance  frequency  of  the  tube  and  then  drilled 
holes  of  various  sizes  in  the  side.  After  each  new  hole  was  drilled,  the  resonance  frequency  was 
checked  and  recorded.  The  result  is  shown  in  table  1 . 
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Figure  11.  Picture  of  the  tube  used  for  varying  the  diameter  of  the  side  branch. 

Table  1.  Result  of  measuring  the  resonant  frequency  of  tubes  with  side  branches  of  various  diameters. 


Diameter  (m) 

(a/b)2 

te(m) 

Lefr(m) 

c/2Leff 

(Hz) 

Measured  Freq. 
(Hz) 

%  Error 

0 

— 

— 

— 

— 

517 

— 

0.00165 

248 

0.004 

1.146 

146 

527 

— 

0.00305 

72.7 

0.00442 

0.475 

352 

553 

— 

0.00561 

21.5 

0.00518 

0.26531 

630 

595 

5.6 

0.0123 

4.47 

0.00718 

0.1861 

899 

745 

17.1 

0.01587 

2.68 

0.00826 

0.17614 

950 

783 

17.6 

0.02222 

1.27 

0.01016 

0.1669 

1003 

842 
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The  most  important  thing  to  notice  is  that  the  diameter  of  the  side  branch  hole  has  a  significant 
effect  on  the  resonant  frequency  of  the  tube.  Also  note  that  when  (a/b)2  >20,  the  error  is  very 
large.  The  reason  for  this  is  that  the  effective  length  of  the  tube  becomes  larger  than  the  actual 
tube,  which  is  0.308  m.  In  terms  of  resonance  in  the  tube,  this  is  not  possible. 

Now,  I  turn  my  attention  to  the  cell  shown  below  in  figure  12.  This  is  an  actual  cell  used  in  a 
photo-  acoustic  spectrometer.  Note  that  this  structure  is  different  than  that  shown  in  figure  10. 
The  microphone  chamber  is,  in  effect,  a  tube  with  varying  diameter.  The  situation  is  shown  in 
figure  13  below.  The  equivalent  matrix  for  this  structure  is  given  in  equation  72. 
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Figure  12.  Diagram  of  actual  acousto-optic  spectrometer. 
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Figure  13.  Structure  of  side  branch  with  two  sections. 


Using  the  dimensions  given  in  figure  13  under  “Microphone  Channels”  in  equation  72  and 
solving  for  the  input  impedance  of  the  side  branch,  we  find  that  the  imaginary  part  of  Zin  is 


lm[Zin]  = 


cos  kLx  sin  kL2-(j^  sin  kLx  cos  kL 


(73) 
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By  inserting  this  term  into  the  bracketed  portion  of  equation  7 1  and  solving,  we  are  led  to  a 
resonance  frequency  of  20,235  Hz.  The  measured  resonance  frequency  was  19,800  Hz,  so  that 
the  error  is  2.1%. 

I  mentioned  earlier  that  there  would  be  a  discussion  of  the  sources  of  error  in  our  calculations. 
One  source  of  error  found  in  the  experiment  with  varying  side  branch  diameters  is  the 
imperfection  of  the  drilled  holes.  The  sharpness  of  the  edges  does  have  an  effect  on  the 
impedance  of  the  hole.  In  our  case,  one  can  see  from  figure  12  that  the  hole  is  very  ragged. 
Another  source  of  error  comes  from  our  assumption  that  the  ends  of  the  tubes  are  an  infinite 
open  with  no  flange  structure.  This  is  not  the  case  in  any  of  the  situations  discussed.  The  final 
source  of  error  to  be  discussed  here  is  the  speed  of  sound  in  the  tube.  In  all  cases,  we  assumed 
perfectly  rigid  and  perfectly  smooth  sides  on  the  tube.  Obviously  this  can’t  be  the  case.  We  could 
and  have  assumed  that  the  formula  given  at  the  end  of  section  2.1  is  valid,  but  in  the  case  of  our 
PVC  tube  this  does  not  seem  to  hold  true.  In  that  case,  the  speed  was  measured  to  be  334  m/s, 
but  this  measurement  was  dependant  on  the  assumption  that  the  end  correction  was  0.3D.  The 
calculated  value  from  equation  16  is  338.76  m/s. 


7.  Conclusions 


We  have  shown  that  it  is  quite  reasonable  to  assume  that  in  our  photo-acoustic  cells  the  acoustic 
waves  travel  as  plane  waves.  This  assumption  is  also  supported  by  Tijdeman  (5),  where  he  shows 
that  the  velocity  profile  of  waves  that  a  have  S  (shunt  wave  number)  greater  than  10  are  fairly  flat 
across  area  of  the  tube.  In  our  case  of  a  1-mm-diameter  representative  cell,  S=44,  which 
corresponds  to  a  velocity  profile  that  is  very  flat  across  the  area  of  the  tube.  We  developed  a 
matrix  method  of  calculating  the  input  impedance  of  the  acoustic  resonators,  and  using  this  can 
predict  the  resonant  frequencies  within  a  few  percent.  This  method  can  be  extended  to  more 
complex  structures  for  side  branches  with  similar  results.  We  also  showed  that  the  diameter  of 
the  side  branch  hole  is  critical  in  the  design  of  a  photo-acoustic  resonant  cell.  For  side  branch 
diameters  that  are  less  than  one-quarter  the  diameter  of  the  main  bore  the  resonant  frequency 
cannot  be  predicted. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


x  equilibrium  coordinate  of  a  thin  sheet  of  molecules 
E,  displacement  of  the  sheet  along  the  axis  of  propagation 
u  velocity  of  the  sheet  u  =  — 

p  instantaneous  density  at  any  point 

po  constant  equilibrium  density  of  the  medium 

5  condensation  at  any  point  defined  as  5  =  -  —  or  p  =  po(l  +s) 

Po 

P  instantaneous  pressure  at  any  point 

Po  constant  equilibrium  pressure 

p  excess  pressure  or  acoustic  pressure  at  any  point  defined  as  p  =  P-Po 
c  velocity  of  wave  propagation 

vp  phase  velocity 

r  propagation  constant  in  Zwikker-Kosten  derivation  of  phase  velocity 
r l  reflection  coefficient  at  the  load  of  a  transmission  line 

SA  thermal  boundary-layer  thickness 

v  kinetic  viscosity 

t]  dynamic  viscosity 

y  ratio  of  the  specific  heats 

Cy  specific  heat  at  constant  volume 
Cp  specific  heat  at  constant  pressure 

k  propagation  constant  everywhere  except  velocity  calculation 

r  radius  of  the  tube 

D  diameter  of  the  tube 

z  propagation  direction 

(o  radian  frequency 
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L  length  in  tube 

Z  impedance 

Y  admittance 

V  voltage 

/  current 

te  effective  length  of  a  side  branch 

Xb  effective  length  of  a  side  branch  corrected  for  the  diameter  of  the  main  tube 
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